clc;
clearvars;
close all;
a  = [1.3565   1.3027   1.2385   1.0476   1.2292   1.1539   1.1691   1.1615   1.0311   1.1466   1.0252   1.0430   1.2112   0.9959   1.0533   1.1535   0.9774   1.1264   1.1036   1.0686   1.1354   1.0377   1.0813   1.1085   1.0541   1.0441   0.9759   0.9958   0.9598   1.0950   0.9332   0.9656   1.2021   0.9859   0.8690   0.9641   1.0369   0.8871   0.8952   0.8617   0.9419   0.9731   0.9381   0.8604   0.8737   0.8209   0.8615   0.8813   1.0104   0.9922   0.8276   0.8299   0.8605   0.8362   0.8599   0.8262   0.8685   0.8676   0.7966   0.8816   0.8489   0.7529   0.8514   0.8050   0.7967   0.9144   0.8381   0.8091   0.9158   0.8165   0.8225   0.8464   0.7913   0.8437   0.8487   0.8451   0.8377   0.8634   0.8066   0.8391];
b = [1.0183   0.9238   0.9184   0.9201   0.9155   1.1478   0.9597   0.9154   1.1238   0.9753   0.9932   1.0375   0.9434   1.0026   0.9969   1.0290   1.0237   0.9010   1.0128   1.0393   0.9490   1.0992   0.9419   0.9273   1.0513   0.8929   0.9288   1.0377   0.8835   0.8743   1.0749   1.0801   0.9005   0.9013   0.9777   0.9461   0.9068   0.8415   0.8999   0.9617   0.9915   0.8965   0.8869   0.8556   0.8425   0.8713   0.9322   0.9753   0.9161   0.8579   0.8272   0.8447   0.8485   0.8123   0.8450   0.7810   0.8421   0.8393   0.8092   0.8582   0.9001   0.8044   0.8540   0.8312   0.7561   0.8337   0.8199   0.7763   0.8288   0.8261   0.8485   0.8310   0.8095   0.8216   0.7967   0.7902   0.8178   0.8454   0.0000];
options.verbosity = 2;
n = numel(a);
options.u_rows = [n];
[U, S, V] = spx.fast.bdhizsqr_svd(a, b, options);
